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ABSTRACT 

The prevailing theory for the origin of cosmic magnetic fields is that they have 
been amplified to their present values by the turbulent dynamo inductive action in the 
protogalactic and galactic medium. Up to now, in calculation of the turbulent dynamo, 
it has been customary to assume that there is no back reaction of the magnetic field on 
the turbulence, as long as the magnetic energy is less than the turbulent kinetic energy. 
This assumption leads to the kinematic dynamo theory. However, the applicability of 
this theory to protogalaxies is rather limited. The reason is that in protogalaxies the 
temperature is very high, and the viscosity is dominated by magnetized ions. As the 
magnetic field strength grows in time, the ion cyclotron time becomes shorter than the 
ion collision time, and the plasma becomes strongly magnetized. As a result, the ion 
viscosity becomes the Braginskii viscosity. Thus, in protogalaxies the back reaction sets 
in much earlier, at field strengths much lower than those which correspond to field- 
turbulence energy equipartition, and the turbulent dynamo becomes what we call the 
magnetized turbulent dynamo. In this paper we lay the theoretical groundwork for 
the magnetized turbulent dynamo. In particular, we predict that the magnetic energy 
growth rate in the magnetized dynamo theory is up to ten time larger than that in the 
kinematic dynamo theory. We also briefly discuss how the Braginskii viscosity can aid 
the development of the inverse cascade of magnetic energy after the energy equipartition 
is reached. 
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1. Introduction 

One of the most important and challenging questions in astrophysics is the origin of strong and 
large-scale magnetic fields in galaxies and protogalaxies. It is now widely accepted that the strong 
cosmic magnetic fields were produced by the turbulent dynamo inductive action driven by the fluid 
motions in a galactic and/or protogalactic medium (Vainshtein & Zel'dovich 1972; Zweibel & Heiles 
1997; Kulsrud 2000). Yet full understanding of all stages of this production has not been achieved. 
There are two alternative theories on how and when the magnetic fields have been produced. The 
first theory, the galactic dynamo theory, also known as the a—Q dynamo theory, states that the 
fields have been primarily amplified in differentially rotating galactic disks after the galaxies had 
been formed (Parker 1971; Vainshtein & Ruzmaikin 1972; Beck et al. 1996). The galactic dynamo 
involves several crucial unsolved problems (Rosner & Deluca 1989; Zweibel & Heiles 1997, Kulsrud 
1999). The main problem is that in the ct-Vl theory it seems to be extremely difficult to expel a 
fraction of the magnetic flux from a galactic disk in order to produce the net magnetic flux (Rafikov 
& Kulsrud 2000). In addition, observations indicate the presence of microgauss magnetic fields in 
galaxy clusters and in early galaxies at high redshifts (Perry 1994; Kronberg 1994). It is hard to 
explain such strong fields by the galactic dynamo theory (Zweibel & Heiles 1997). In this paper we 
accept the second theory for the origin of cosmic magnetic fields, the primordial dynamo theory, 
which states that the galactic and extragalactic magnetic fields have primarily been produced in 
protogalaxies, i. e. before the galaxies were formed (Pudritz & Silk 1989; Kulsrud & Anderson 
1992; Kulsrud et al. 1997; Kulsrud 2000). Of course, these fields were subsequently modified in the 
rotating galactic disks after the galaxies were formed. 

In order to understand how the magnetic fields can be built up in protogalaxies, let us briefly 
discuss the physical conditions that were present there 1 . Let us assume the following typical 
values for the total mass M of a protogalaxy, the total to baryon mass ratio £, and the protogalaxy 
size L: M ~ 1O 12 M , £ ~ 10, and L ~ 0.2 Mpc. Then, the number density of the gas in the 
protogalaxy is n ~ £~ 1 ML _3 m~ 1 ~ 5 x 10~ 4 cm -3 oc £ -1 ML -3 (m p is the proton mass, for 
convenience, in this paragraph we give the scaling of physical parameters with the three "basic" 
parameters £, M and L). Assuming the energy virial equilibrium in the protogalaxy, we easily 
estimate the gas temperature, T ~ k^Gm v ML~ x ~ 2 x 10 6 K oc ML^ 1 (G and fee are the 
gravitational and the Boltzmann constants). This temperature is very high, while the density is 
very low. As a result, the gas is fully ionized, and the viscosity is dominated by ions, not by 
neutrals. The ion collision time is t{ ~ 20T 3 / 2 n~ 1 A~ 1 sec ~ 3 x 10 12 sec oc ^M 1 / 2 ^/ 2 (in this 
formula the temperature is in degrees K, the density is in cm" 3 , and A c ~ 30 is the Coulomb 
logarithm assumed to be independent of £, M and L, Braginskii 1965). The virial thermal speed 
is Vt ~ (2^T/m p ) 1 / 2 ~ 2 x 10 7 cm/s oc M 1 / 2 L" 1 / 2 . The ion kinematic viscosity can be estimated 
as v = 0.96 k B Tti/m p ~ 5 x 10 26 cm 2 /s oc £M 3 / 2 L 1 / 2 (Braginskii 1965). The Spitzer resistivity is 
rj s = 6.53 x 10 12 T~ 3 / 2 A C cm 2 /s ~ 8 x 10 4 cm 2 /s oc M~ 3 / 2 L 3 / 2 (in this formula the temperature is 



1 The numerical quantities given below refer to common values at the time (redshift) when the protogalaxies form. 
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in degrees K, and the Coulomb logarithm is assumed to be constant, A c ~ 30, Spitzer 1962). Now 
we estimate the Reynolds and the Prandtl numbers, R ~ Vr/kov ~ 10 4 oc ^ _1 M _1 (ko = 2ir/L 
is the minimal wave number in the protogalaxy) and Pr ~ v/rj s ~ 10 22 oc ^M 3 L _1 , they are 
very large. The viscous cutoff scale of the turbulence can be estimated as 2irk~ 1 ~ i?~ 3 / 4 L ~ 
10 _3 L oc £ 3/4 M 3 / 4 L, the resistive cutoff scale for the magnetic field as 2itk~} ~ Pr^^R'^L ~ 
1(T 14 L oc ^/ 4 M" 3 / 4 L 3 / 2 , and the ion mean free path as A; ~ R^L ~ W~ 4 L oc (ML. Thus, 
there is a hierarchy of scales in the protogalaxy, L S> 2irk~ 1 S> A« 3> 2irk~\ Therefore, we can 
use the single-fluid magnetohydrodynamic (MHD) equations for the description of plasma, and can 
consider the plasma to be nonresistive and incompressible on scales L > 2i:k~ l > \. 2 

It is very important that in a protogalaxy the ion cyclotron period in the magnetic field, 
a;" 1 = m p c/eB ~ 10~ 4 -B -1 G • sec, is shorter than the ion collision time, tj ~ 3 x 10 12 sec, provided 
that the field strength is larger than ~ 10 -17 -10~ 16 G. On the other hand, the magnetic energy 
becomes comparable to the kinetic energy of the smallest turbulent eddies (which are on the viscous 
cutoff scale) if the field strength exceeds ~ (47rm p nV^i? -1 / 2 ) 1 / 2 ~ 10~ 7 G. Thus, there is a 
very broad range of magnetic field strengths, at which the magnetic pressure and tension are still 
negligible, while the presence of the field is already important. This is because the plasma is 
strongly magnetized, uj^ 1 <ti ti, and the magnetic field controls the microscopic motions of ions, so 
that the viscous forces are given by the Braginskii viscous stress [Braginskii 1965; see also eq. (1)] 
and are different from the viscous forces in a field-free plasma. 



2 Note that the Debye length (fcsT /2-ime 2 ) 1 ^ 2 ~ 6 x 10 5 cm oc £ 1//2 L is extremely small in protogalaxies. Also we 
can consider the plasma to be incompressible on scales < L because the plasma velocities at the largest scales, ~ L, 
are of the order of the sound speed, and all velocities at smaller scales are smaller. For estimation purposes, we can 
use the MHD equations even for scales <C A;, if we reduce the molecular viscosity by factor (fcA;) -1 <C 1. 
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Fig. 1. — There are five major stages of the production of galactic and extragalactic magnetic fields, 
as the field strength grows up from zero to < 1CT 5 gauss. In this paper we consider the magnetized 
turbulent dynamo stage in protogalaxies, shaded in the plot. 
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We believe that there are five major stages of the production of the strong cosmic magnetic 
fields, see Figure 1. During the first stage, in a protogalaxy undergoing gravitational collapse, the 
Biermann battery builds a seed magnetic field linearly in time, on a time scale approximately equal 
to the free-fall time, ~ 1 billion years. The resulting seed field is of the order of ~ 10~ 18 G on the 
viscous scale (Pudritz & Silk 1989; Kulsrud et al. 1997; Davies & Widrow 2000). 

During the second stage, when the plasma is unmagnetized, uj^ 1 S> U, the seed field is expo- 
nentially amplified by the kinematic turbulent dynamo inductive action. The kinematic dynamo 
builds the field up to approximately 10~ 16 G, when the plasma becomes magnetized, uj~ x <C t«. The 
time scale of the kinematic dynamo is very short, it is of the order of the smallest eddy turnover 
time (Kulsrud et al. 1997), ~ 10 million years, which is much smaller than the gravitational col- 
lapse time of the protogalaxy, several billion years. The main assumption of the kinematic dynamo 
theory and, basically, the strict definition of it, are that the growing magnetic field stays so weak, 
that it does not affect the fluid motions, i. e. that there is no the back reaction of the field on the 
turbulence. The results we obtain in this paper reduce to those of the kinematic dynamo theory in 
the limit of very weak field strengths (i. e. when the plasma is unmagnetized, and there is no the 
back reaction). 

The third stage starts when the field grows above ~ 10~ 16 G, the ion cyclotron time in 
the magnetic field becomes shorter than the ion collision time, and the plasma becomes strongly 
magnetized. During this stage the magnetic field strongly affects the dynamics of the turbulent 
motions on the viscous scales 3 , by completely changing the viscosity (see Section 2), despite the 
fact that the magnetic energy is still very small compared to the kinetic energy of the turbulence! 
We call this stage as the magnetized turbulent dynamo. In previous theories this stage has not been 
recognized. The main goal of this paper is to construct a theoretical model for it. The time scale 
of the magnetized dynamo is about the same as the kinematic dynamo time scale, ~ 10 million 
years. 

So far the field scale is of the order of the viscous scale or less, and the magnetic field is 
incoherent in space. The fourth stage starts when the magnetic field strength grows above ~ 
10~ 7 G. The field energy becomes comparable to the kinetic energy of the smallest turbulent 
eddies, and the Lorentz forces become dynamically important in the plasma. During this stage 
the turbulent motions are dissipated by the growing field, the turbulent energy spectrum becomes 
truncated at larger and larger scales, and the turbulent energy is eventually transferred into a large- 
scale strong magnetic field with its energy comparable to the kinetic energy of the fluid motions 
on the largest scales in the protogalaxy (which is approximately the same as the thermal energy). 
This process (which is still under debate) is called the inverse cascade (Vainshtein 1982; Kulsrud 
& Anderson 1992; Beck et al. 1996; Kulsrud 2000). We discuss it qualitatively in Section 5. The 



3 The turbulent motions on the viscous scales are the most important in the dynamo theory. This is because the 
magnetic field is primarily amplified by the turbulent eddies on these scales. These eddies are the smallest ones, they 
have the shortest turnover times and produce the largest velocity shearing. 
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time scale available for the inverse cascade process is of the order of the largest eddy turnover time, 
~ 1 billion years. 

The turbulent dynamos and the inverse cascade may not have time to amplify the field up to 
microgauss values, which are observed in galaxies. A crucial question is how far they go? This 
paper addresses this question and is concerned with the rate of the magnetic field built up by the 
magnetized turbulent dynamo. 

Finally, the fifth stage is the galactic dynamo, which happens in the differentially rotating 
galactic disc after the galaxy is formed. This process modifies the strong field that was initially built 
up in the protogalaxy, on a time scale of the order of the rotation time of the galaxy, ~ 300 million 
years. The galactic dynamo theory is not discussed in this paper. 

In Section 2 we formulate the basic equations of the magnetized dynamo theory. In Section 3 we 
calculate the statistical correlation functions for turbulent velocities, V, in a strongly magnetized 
plasma by making use of the quasilinear expansion in time of the MHD equations for both the 
velocities and the magnetic field, similar to the expansion used by Kulsrud and Anderson (1992). 
We find that, contrary to the Kolmogorov velocities, the turbulent velocities in the magnetized 
plasma are strongly anisotropic on the viscous scale, as one might expect, because the magnetic 
field sets "a preferred axis in space" . In our calculations we make two working hypotheses. First, we 
assume that for the purpose of magnetic energy calculation, the tensor b a p = b a bp, where b is the 
field unit vector, can be taken to be constant in space in the beginning of the expansion in time. This 
is our first hypothesis, which basically relies on our assumption that in the magnetized turbulent 
dynamo case the magnetic field has a folding structure similar to the one that exists in the kinematic 
turbulent dynamo case (see Figure 2; Maron & Cowley 2001; Schekochihin et al. 2002). Second, 
we find that there are velocity modes which are not damped by the Braginskii viscous forces and, 
therefore, are divergent unless we incorporate the MHD non-linear inertial terms into our quasilinear 
expansion. We assume that these non-linear terms, which limit the divergent velocity modes by 
coupling them to other modes, may be included into our theory by allowing for rotation of velocity 
vectors relative to the magnetic field unit vectors 4 . This is our second hypothesis. In Section 4 
we use the correlation functions for the turbulent velocities found in Section 3 to calculate the 
evolution of the magnetic energy in the magnetized turbulent plasma. We start with calculations 
of the total magnetic energy growth rate in Subsection 4.1. In Subsection 4.2 we derive the integro- 
differential mode coupling equation for the magnetic energy spectrum. Our mode coupling equation 
is a more general version of the corresponding equation of Kulsrud and Anderson (1992), obtained 
for the kinematic dynamo. In Subsection 4.3 we consider the magnetic energy spectrum on small 
subviscous scales. On these scales the mode coupling equation greatly simplifies and becomes a 
homogeneous differential equation. Finally, in Section 5 we give our conclusions. We also discuss 
the peculiarities of the inverse cascade in a strongly magnetized turbulent plasma. 



4 This rotation is essentially due to Coriolis forces that make the velocity rotate differently than the field direction. 
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2. Basic Magnetized Dynamo Equations 

Hereafter we consider the magnetized turbulent dynamo stage, the shaded region in Figure 1. 
During this stage the plasma is strongly magnetized, Wjtj S> 1, but the magnetic energy is still 
small compared to the turbulent kinetic energy, and therefore, the magnetic Lorentz forces can be 
neglected. The viscous forces acting on turbulent velocities V in a strongly magnetized incom- 
pressible fully ionized plasma are determined by the Braginskii viscosity stress tensor (Braginskii 
1965), 

= -y{^b a bp - Sa^bJtvV^ , (1) 

where b = ~B/B is the unit vector along the magnetic field. Note, that this stress tensor depends on 
the field unit vector, but is independent of the magnetic strength (as long as the plasma is strongly 
magnetized, and turfi 3> 1). Thus, during the magnetized turbulent dynamo stage the magnetic 
field strongly affects the turbulent motions on the viscous scales by changing the viscous forces, 
even though the Lorentz forces are negligible. 

The MHD equations for the turbulent velocities V in an incompressible strongly magnetized 
plasma, neglecting the Lorentz forces, are (Landau & Lifshitz 1984) 

d t v a = -p[ a + f a -^ p -{v a v p ) t p 

= -P!L + fa + MkbMuV^),fs-(V a V p ), p , (2) 

V a , a = 0, (3) 

where f is the force driving the turbulence, and P' is the hydrodynamic pressure. Here and below we 
always assume summation over repeated indices. In order to shorten notations, we use d t = d/dt, 
and spatial derivatives are assumed to be taken with respect to all indices that are listed after " , " 
signs 5 . To obtain the second line of equation (2), we use formula (1) for the viscous stress and 
incorporate the isotropic part of the stress into the pressure P" . 

It is difficult to solve equations (2) and (3) directly because they are very complicated. We 
also do not know the exact expression of the driving force f , but its statistics is the same as for 
an unmagnetized plasma. Therefore, let us proceed as follows. We represent f by introducing 
subsidiary incompressible turbulent velocities U which, by definition, satisfy equations 

d t U a = -P^ + f a + (l/5)uAU a -(U a U^)^, (4) 
U a , a = 0, (5) 

where AU a = U a $p . (U will be essentially the Kolmogorov turbulent velocities.) 

Let us analyze and compare equations (2) and (4). First, if for the moment we formally average 
the Braginskii viscosity term 2>v{b a bpb^b y V^^) ^ in equation (2) over all directions of an isotropic 



5 For example, (V a Vg), 7 = Vp(dV a /dx~/) + V a (dVp/dxj), and V Qi) 9 7 = d 2 V a /dxpdx~t. 
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magnetic field, then it reduces to (l/5)vAV a , which coincides with the isotropic viscosity term in 
equation (4). Therefore, v e s = [l/h)u could be considered as an effective reduced viscosity for an 
incompressible fully ionized plasma in the presence of a magnetic field that is isotropically tangled 
on subviscous scales. In other words, the Braginskii viscous forces "are doing worse" at dissipating 
the turbulent motions, as compared with the standard isotropic viscous forces in a field- free plasma. 
Second, note that equations (2) and (4) have the same driving force f . By taking the driving force 
to be the same, we assume that this force comes from larger turbulent eddies. These larger eddies 
are on scales larger than the viscous scales, and therefore, these eddies "do not know" whether the 
viscous forces are of the Braginskii type or of the standard isotropic type. 

Now, note that equation (4) is a familiar hydrodynamic equation with a standard isotropic 
viscosity term. However, is has a reduced molecular viscosity, (l/5)u instead of v. Therefore, we 
assume that the solution of equations (4) and (5) is the incompressible, homogeneous, isotropic and 
stationary Kolmogorov turbulence with the effective reduced viscosity 

u eG = (l/5)i/. (6) 

As a result, the statistics of the Fourier coefficients of the turbulent velocities U, 

1 f L ' 2 

-L/2 



U lta (u) = -= UteWe^dt, £/ kQ ( t )= / U a (t,r)e-* kr d 3 r, (7) 

\/ 2ir J-oo L J —L/2 



is given by the following formulas 6 : 

(CM")) = 0, (8) 
(Uka(u)Uvp(u')) = (Ul^ a {-u)Uup{J)) = Ju k 8^8u^5(J + u), (9) 

S af3 = f <W - k a k(3, (10) 

Ju,k = J ofc (l + rV) _1 , (11) 

r(fc) = r(0) (k/k )- 2/3 = (l/k Uo) (k/k y 2/3 , (12) 

j ~ J (Uo/^k )(k/k y 13 / 3 , k <k<k u = (5U /k v) 3 / 4 k , (u) 

° k ~ 1 0, k<k = 2vr/L, k > k v . 

Here and below (...) means ensemble average over all realizations of the turbulence, 5(to' + lo) 
is the Dirac <5-function, 8 a n and <5k',-k are the one- and three-dimensional Kronecker symbols 
respectively, k = k//c is the unit vector along k, Uq ~ Vt is the largest eddy velocity, k^ = 2ir/L 
is the smallest wave number of the turbulence, and k v ~ (Uo/kou e g) 3 / 4 ko = (bUo/kov) 3 / 4 ko is the 
viscous cutoff wave number of the turbulence 7 . In equation (9) we keep only the normal part of 



6 Note that these formulas are similar to those of Kulsrud et al. 1997, but differ by numerical coefficients from 
those of Kulsrud & Anderson 1992. 

7 k l , is determined by the balance between the inertial and the viscous terms of eq. (4), l/r(k v ) ~ v e sk^,. 
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the turbulence and drop the helical part, since the latter is negligible on the scales of the smallest 
turbulent eddies, which are the principal drivers of the field evolution (Kulsrud & Anderson 1992). 
To obtain equation (11), we assume that the time correlation function of the turbulent velocities 
has an exponential profile 8 , 



<tWt)Ck/j(0> = ^ _M ' l/T ^k',-k, (14) 



where r is the eddy decorrelation time given by equation (12) for the Kolmogorov turbulence. Note 
that the averaged total kinetic energy of the fluid motions, per unit mass, is 



\(m, r)] 2 ) = \ £<|U k (*)| 2 > = lj2^jL = l£ I(k) dk = \ul (15) 



where the Kolmogorov energy spectrum is I(k) = 4irk 2 (L/2ir) 3 Jok/T = (2/3) (C/q /k )(k/k ) 5 / 3 if 
k € [ko, k u ], as it should be (Kulsrud et al. 1997). 

Now let us subtract equations (4) and (5) from equations (2) and (3) to eliminate the unknown 
driving force f , and let us introduce the back-reaction velocity v d = V — U. We have 

d t v a = -P, a + ^v{bap^v^ v + b a p^U^ v )^-{\/h)v/\U a -{v a Up + U a Vi3 + v a Vi3)^, (16) 

V a ,a = 0, (17) 

where the pressure P = P" — P'" . Here and below we use the following symmetric tensors: 

&q/3 7 5 = babpb-yh , b a/3y = f bj)^ , b af3 d = bab/3 . (18) 

Velocity v, which satisfies equations (16) and (17), can be considered as the correction to the 
Kolmogorov velocity U. This correction is non-zero only on the viscous scales, and results from 
the strong back reaction of the field on the turbulence via the Braginskii viscosity tensor (1). 

Finally, the MHD equation for the magnetic field B is (Landau & Lifshitz 1984) 

d t B a = V^Bp-VpB^p, (19) 

where the plasma velocities V = \J-\-v are incompressible, and we neglect resistivity. Consequently, 
the equations for the magnetic field squared, B 2 , and for the magnetic field unit vector, b, are 

8 t B 2 = 2V a ^B a Bp - Vp(B 2 )^ , (20) 
d t K = V a ^bp - Vs i7 6 a/ 3 7 - Vpb^p . (21) 



8 Using a Gaussian time correlation profile, (U(t)U(i')) tx e~'*~* ' / 2t , would be more appropriate. In this case 
equation (11) would become J^ k = Joke~ T w ^ 2 ■ However, we prefer the exponential profile because it is easier to 
deal with (e.g., for Gaussian integrals it is not possible to close integration contours at infinity in the complex plane). 
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3. Statistics of Turbulent Velocities in Strongly Magnetized Plasmas 

In order to find the evolution of the magnetic field B, we must derive the correlation functions 
for the total velocities V, resulting from the Braginskii viscosity. We calculate these velocity 
correlation functions in this section. 

Let us assume that we know the magnetic field at zero time, B| t=0 = °B(r) and b| t=0 = °b(r), 
and that the back-reaction velocity v is initially zero, v\ t=0 = °v(r) = 0. 9 Then we advance 
the magnetic field and the back-reaction velocity to some future time, t > 0, by the nonlinear 
terms, i. e. by integrating equations (16), (19) and (21) twice in time. This quasi-linear expansion 
procedure is similar to the calculations of Kulsrud and Anderson (1992). Considering t as the 
expansion parameter 10 , up to the second order, we have 



B(t,r) = °B(r) + 1 B(t,r)+ 2 B(t,r), (22) 

b(t,r) = b(r) + 1 b(t,r) + 2 b(t,r), (23) 

v(t,r) = 1 v(t,r)+ 2 ^(t,r), (24) 

V(t, r) = 1 V(t, r) + 2 V(t, r) = [U(t, r) + 1 v(t, r)] + 2 v{t, r). (25) 



Here, V = XJ+v is the total fluid velocity, and the Kolmogorov turbulent velocities U are considered 
to be given and to be of the first order (Vainshtein 1970, Kulsrud & Anderson 1992). 

Now, we substitute the above expansion formulas into equations (16)-(21). We find that the 
zero order equations are 



d t °B a = 0, d t °b a = 0, V = 0, °V a = 0, (26) 
the first order equations are 

d^Ba = ^aA-V^. (27) 

dt 1 b a = 1 y r a,/3°b/3 — 1 Vp^ b a/ 3 y — 1 y / 3°6 Q , ij a , (28) 

dt x v a = - 1 P,a + M K^v 1 Vw),p + M K^vU^) !/3 -(l/5)vAU a , (29) 

= 0, (30) 

l V a = U a + l v a , (31) 

and the second order equations are 

dt 2 B a = 1 V aj p 1 B /3 + 2 V at p Q Bp - l V p l B a> p - 2 V p°B a ^ , (32) 

dt 2 b a = 1 Va )l 3 1 b l 3 + 2 V ayj 3°bp - 1 V / 3^ 1 b a p 1 - 2 V "/3 )7 °6 a/ 3 7 - 1 V / 3 1 b a ,/3 - 2 V(3°b a)l 3 , (33) 



9 A nonzero initial back-reaction velocity would lead to transients, which would be dissipated anyway. 

10 To be more formal, we need to introduce a dimensionless variable £ = t/At, and to consider U a> gAt and V a ,pAt, 
which are dimensionless, as the expansion parameters. 
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dt 2 V a = - 2 P, a + 3u[°b a ^ u 2 V^ + ^^V^u + - + + (34) 

^«,a = 0, (35) 
2 V a = 2 v a . (36) 

Here, of course, the pressure is also expanded, P = l P + 2 P. 

First, let us solve the first order equations (29) and (30). We Fourier transform these equations 
in space and in time, r — > k and t — > oj, by making use of the discrete and the continuous Fourier 
transformations respectively [see eq. (7)]. We have 

(-iuj + Q rd ) 1 v ka (u;) = -ikJPk + 3vikp ^ ik^v^^ + U V/l ] °6 k // a/3/li , + (l/5)vk 2 U ka , (37) 



k' 

k"=k-k' 



fca^kaCw) = 0, (38) 

where 1 v^(uj), 1 Pk(^), Uk(w) and bk a /3nv are the Fourier coefficients. In the left-hand-side of 
equation (37) we add a damping term, Q^v^, in order to account for the rotation of velocity 
vectors relative to the magnetic field unit vectors, see the discussion on page 15. We estimate the 
effective rotational damping rate O r d on page 17, and show it is a constant. In general, it may be 
a function of k = |k| and of (b • k) 2 . 11 

Now, we multiply equation (37) on the left by tensor S^ a = Sj a — kjk a to eliminate the pressure 
term by making use of the incompressibility condition (38). Interchanging indices, and using the 
symmetry of tensor °bk"a/3ni> with respect to its spatial indices, we obtain 

M\a a yp 1 V\ !i ip = iFka , (39) 

°M ka yp = {-iio + ^ rd )(5 kik /(5 a/ 3 + Si/S^kuk'^bw^p , k" = k-k', (40) 

^ka = [-°-A/fka,k'/3 + (~tW + n Td + uk 2 /5) 5 k ,k'<W] , (41) 

where we use convenient matrix notation, so that summation is implicitly assumed over repeated 
spatial indices and wave numbers. The matrix operator Q M. ka ^ p(uj) is of the zero order, while 
"the driving force" 1 ^ r ka(<^') is of the first order. If there exist inverse matrix °.M kc | k /«, then, using 
equations (31), (39) and (41), we obtain the Fourier coefficient of the first order total velocity 1 V, 

l V ka {io) = U ka + l v< ka = (-ito + n Td + vk 2 /b) °7Wkik'/3 U V p(w) . (42) 



n It depends on the square of b • k because the Braginskii viscosity is invariant under reflection b — > b (gyrating 
ions "do not care" about the exact direction of b). 
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Fig. 2. — The folding structure of magnetic fields (for simplicity shown in two dimensions). The 
field is nearly straight and strong in Region I. The field is very curved but weak in Region II. 
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Now, let us consider an important case for which we can invert matrix vWka,k'/3) given by 
equation (40). Let (°B-V)°B = 0, i. e. let at zero time the magnetic field vary only in the direction 
perpendicular to itself. This is equivalent to the magnetic field lines being initially straight (no 
curvature), and ®b a p = °6 Q °6 j a = const. This model of straight magnetic field lines is not as 
artificial as it seems at first glance, because of the following arguments. First, it is known that 
in the kinematic turbulent dynamo case the magnetic field lines have a folding pattern, shown in 
Figure 2, (Schekochihin et al. 2002; Maron & Cowley 2001). This folding pattern implies that in 
the bulk of the volume the field is strong and has small curvature, k± 3> k\\ ~ k v (fcy and k± are the 
wave numbers parallel and perpendicular to the field lines), while in a small fraction of the volume 
the field is weak and curved, k± ~ fey 3> k v . The regions of weak and curved field, Region II in 
Figure 2, can be disregarded as long as we consider the volume averaged magnetic field energy and 
are not interested in the field curvature. As for the regions of strong magnetic field with small 
curvature, Region I in Figure 2, the field lines in these regions can be well approximated by the 
straight field lines on scales k which satisfy k > kit ~ k v . Now, we make assumption that even for 
the magnetized turbulent dynamo case the magnetic field has a folding structure similar to that 
for the kinematic turbulent dynamo case. This assumption will be our first working hypothesis. 
It is based on the simulations of Maron and Cowley (2001), who found good indications for the 
magnetic field folding pattern in their numerical simulations of MHD turbulence with the Braginskii 
viscosity. Unfortunately, our calculational methods are not adequate to theoretically justify our 
first working hypothesis because of complications of the field curvature calculations. However, there 
exist a simple reductio ad absurdum theoretical argument, supporting the hypothesis, which is as 
follows. If the magnetic field were not folded, i. e. if k± ~ fey > k v everywhere in space, then the 
field would be isotropically tangled on viscous and subviscous scales. In this case the Braginskii 
viscosity could be averaged over the field direction, and it should quickly reduce to the isotropic 
effective viscosity v c r = v/5 [see eq. (6)]. As a result, in this case the magnetized dynamo should 
possess the properties of the kinematic dynamo with this effective viscosity, and would develop the 
folding structure of the magnetic field lines. 

Note that at a given fixed spatial point of the Eulerian system of coordinates the field curvature 
changes in time. Thus, even if the field is initially straight at this point, it may not remain straight 
in the future. In other words, as time goes on, sometimes the spatial point belongs to a region 
I and sometimes to a region II (this reflects the time intermittency and spatial convection of the 
curvature). However, we are not interested in the magnetic energy evolution at each point of space, 
but instead, we are interested in the evolution of the total (spatially averaged) magnetic energy. 
Therefore, at any arbitrary moment of time we consider only those points of space that belong 
to all regions I at this time moment. These points give the dominant contribution to the total 
magnetic energy, and at these points the field lines can be considered to be straight. We apply 
our mathematical model of straight field lines to these points, and we find the change (the time 
derivative) of the total magnetic energy at the time moment we consider (see Section 4) . Because 
the moment of time is arbitrary, our final results give the correct evolution of the total magnetic 
energy in time. 
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Thus, from our first working hypothesis °b a p = const, we have °6k"7^/3 = Sw 1 fl b^b^b^bp, 
and can easily invert matrix (40), 



-l 

ka,k'/3 



3k,k' 



-iuj + Q r d 



$a/3 - 



-iuj + r2 r( j + 3vk 2 n 2 (l — n 2 ) 1 — fi 2 



def oc f 
|U = b • k. 



We substitute this formula into equation (42) and obtain 



— + Q ld 



! -// . . -iuj + Q ld + n 5^ °b 7 °bp ~ 

= , , , n , 00 i ..2 ^(w) 



-iw + fi rd + 20 1 - ,u 2 
where, we introduce the following notations for the viscous damping frequencies 

Cl = (l/5)uk 2 = u cS k 2 , n = (3/2)i/fcV(l - fi 2 ) = (15/2)fi/i 2 (l - / u 2 ). 



(43) 
(44) 

(45) 
(46) 

(47) 
(48) 



The frequency fi depends on because of the anisotropy of the Braginskii viscous stress tensor. 
The frequency Cl is equal to Q averaged over fj,, and represents the averaged rate of the Braginskii 
viscous dissipation [see eq. (6)]. 

Next, we calculate < 1 V" ka (t)) and C^ka W^k'MO), which 

are the ensemble averages of the 
1 V's over all possible realizations of the turbulent motions. First, using formulas (8), (9) and (45)- 
(47), we find 



(^kaH) = 0, 



(49) 



H F (u; O, + tt vd , fi rd ) Sat 



s± x± Ou Ot ' 
°a~i p T °1 °T 



H F (u;q 1 ,q 2 ) = 



+ h f (uj; n + n Td , 2n + n rd 



drf LO 2 + q\ 



s± s:± Ot Ot 
a-y l3r °1 °t 

1-fi 2 



2 ' 



(50) 
(51) 



Note that ( 1 ^ka( LJ ) 1 ^k'/3( u;/ )) = 0) this implies that 1 V k and 1 V k are perpendicular on average. 
Second, we apply the inverse Fourier transformations in time, u,u' — > t,t', to equations (49), (50). 
We have for the required correlation functions of the total turbulent velocities 



(52) 
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= & 



Hp(t — t ; + O rd , O rd 




+ H F {t - *' ; o + o rd , 20 + o rd 7 r T I — 

1 — 



(53) 



where the function 



H F (t ~ f ; 51, (B) = ^ T ^H"; 91, ») e-^-^cL, = ^ /°° 4^4 ^^fP 1 (54) 
2vry_ 00 2vr7_ OCJ u; 2 +^ 1 + r 2 ^ 2 

is the inverse Fourier transform of function J^Hp^u; qi, (72), and depends only on the absolute 
value of the time difference t — t'. In equation (54) we use formulas (11) and (51). 

The main unknown quantity is the effective rotational damping coefficient O rd , which we 
introduced in equation (37) as the O^ 1 ^ term. Let us consider it in more detail. We start with 
the calculation of the ensemble average of 1 V squared 12 . We have 

(V) = (V a (t,r)V a (t,r))=^( 1 y ka (t) 1 F k , a (t))e i ( k + k > = ^(| 1 V k | 2 ) (55) 

k,k' k 

[see eq. (53)]. Using equations (53) and (54), we obtain 

(rv k | 2 ) = tf F (o ; o + o rd ,o rd ) + tf F (o ; o + o rd ,20 + o rd ) 

Jpk Jpk (O + Ord) 2 - Q 2 d Jpk (0 + Ord) 2 ~ (20 + O r d) 2 

r 2 O rd (l + rO rd ) 2 (20 + O rd )[l + r(20 + O rd )] ' 1 ' 

According to equations (55) and (56), we see that if we set the rotational damping rate O rd to zero, 
then the ensemble average of the first order velocity squared would become infinite, ("*" V 2 ) — ► oo as 
O r d —> ! Let us try to understand this divergence problem, and see how we can avoid it. 

First, following the derivation of formula (56) from equation (53), it is easy to see that the diver- 
gence of Cv 2 } occurs only due to the divergence of the 1 V ka modes. In other words, (| 1 V k | 2 ) — ► oo 
as O rd — > 0, while (|v k | 2 ) stays finite 13 . Second, let us refer to equations (46) and (47). On 
one hand, we have &a ^ ka = and k^V^ = 0, as it should be because plasma velocities are 
incompressible. On the other hand, we have °&a ^ ka = and & Q ^ ka ^ 0. Thus, the divergent 
velocity modes, 1 V k , are perpendicular to both vectors k and °b. At the same time, the other, 
non-divergent modes, 1 V k , have nonzero components along °b. Third, let us calculate the ensemble 



12 Note that ( 1 V) = because there is no preferred direction. Of course, there is a preferred axis in space, which 
is along the magnetic field unit vector, see also the footnote 11 on page 11. 

13 Of course, in the degenerate cases, when k _L °b (i. e. fi 2 = 0) or k || °b (i. e. 1 — fi 2 = 0), both 1 V] C and x Vk 
modes become infinite as fLd — » 0. 
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averaged Braginskii viscous dissipation into heat (Braginskii 1965). We have 

k,k' 

- -3z.^A;V ^( 1 ^L 1 ^-k 7 ) ^= -2^^^F(0;O + ^ rd ,2O + O rd ), (57) 

k k 

where we keep only the first order velocities, and make use of formulas °b a p~ /T = const, b a 1 v' ka = 0, 
and of equations (1), (44), (48), (53). We see that only the 1 V k velocity modes are dissipated by 
the Braginskii viscous forces, and the dissipation is proportional to 20 = 3iyk 2 fi 2 (l — n 2 ). 

As a result, we conclude that the divergence of the first order velocities happens because the 
1 V k velocity modes, which are perpendicular to °b, are not damped by the Braginskii viscous 
dissipation. On the other hand, it is clear that the 1 V k modes can not be infinite. What are the 
physical mechanisms which limit them? To answer this question, let us note that velocity modes 
are non-linearly coupled with each other via the inertial term, (V • V)V, of the MHD equation (2). 
This non-linear coupling transforms the divergent 1 V k velocity modes into other modes, which are 
then viscously dissipated. This transformation can be viewed as a continuous rotation of velocity 
vectors relative to the magnetic field unit vector b. Indeed, at a given point of space we can go to 
a reference frame that rotates together with b. In this rotating frame there exist Coriolis forces, 
which act on velocity vectors and force them to rotate relative to the non-rotating field vector. 
These Coriolis forces are caused by the non-linear coupling of velocity modes via the inertial term. 
As a result, a divergent velocity mode, perpendicular to b and not viscously damped, eventually 
rotates out of its initial direction and is transformed into a damped mode. We call this process 
the "effective rotational damping". Of course, it operates only on the viscous scales, on which the 
Braginskii viscous dissipation is significant. On larger scales the viscous dissipation is small and 
the rotation of velocities does not make any difference. 

We could regard f2 rd as a parameter to be determined by numerical simulations. Nevertheless, 
it is of interest to attempt to estimate it from a physical argument. First, the square of the angular 
velocity of the rotation can be estimated as uj 2 ot ~ (l/3)r~ 2 , where t„ is the velocity decorrelation 
time on the viscous scale. The factor 1/3 in this equation comes from the fact that only one of the 
three angular velocity components contributes to the deviation of a divergent velocity mode from 
its original direction perpendicular to the field unit vector b 14 . Second, the typical Braginskii 
viscous damping rate is clearly about \jr v ~ v c ^k 2 , and is larger than the angular velocity of 
the rotation. Let assume it is much larger, uj 2 ot T 2 ~ 1/3 <C 1. Then we can suppose that after 
the divergent velocity mode, in a time interval At, rotates by an angle A(f> relative to its original 
direction perpendicular to b, only the projection of the velocity mode on the plane perpendicular 
to b survives, and all other velocity components are immediately viscously dissipated. As a result, 



This is the component along the vector product b x V. 
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the effective rotational damping rate, fi r d, can be estimated as follows. 

dV AV F(cosA0-l) V Ajl2 V . , 2 At 

~& ~ AF At 2A^ A</> ~ " 2A^ ^ Tv) ~ v 

~ -(l/6)r i 7 1 y~(l/6)i/ cff fc2y~(l/30)^y = _^ rd y 5 ( 58 ) 

ft rd = (l/6)i/cfffe2 = (1/30)^. (59) 

To obtain the last result in the first line of equation (58), we use the random- walk approximation 
for the estimate of Ac/) 2 . 15 

Let us summarize our discussion of the effective rotational damping. As we said, this physical 
damping is associated with the non-linear coupling of velocity modes, which leads to the rotation of 
velocities relative to the magnetic field vectors. The effective rotational damping is very important 
because it limits the velocity modes which are perpendicular to the magnetic field vectors, and 
therefore, are undamped by the Braginskii viscous forces. The non-linear coupling of the veloc- 
ity modes is hard to deal with directly. In particular, the MHD non-linear inertial terms do not 
appear in our first-order equation (37). As a result, to avoid the divergence problem for the un- 
damped velocity modes, we incorporate the non-linear mode coupling and the associated rotational 
damping into our equations in a simple way, as the fi r d 1 ^k damping term in the left-hand-side of 
equation (37). This is our second working hypothesis. Note, that the r2 r d 1 t)k term is isotropic, and 
therefore, it damps not only the divergent velocity modes, perpendicular to the field vector, but all 
velocity modes. This is not a serious problem though, because the rotational damping is smaller 
than the Braginskii damping by a factor ~ 1/6 [see eq. (59)], and our results should be valid 
within a factor of order unity. Also note that on scales larger than the viscous scale the rotational 
damping does not operate. However, on those large scales the turbulence is Kolmogorov, V = U, 
and the back reaction velocities v are zero anyway 16 . 

So far we have considered only the first order velocity, X V, and have found its statistics, given 
by equations (52)— (54) and (59). In order to find the second order velocity, 2 V, we need to solve 
the complicated second order equations (34)-(36). Fortunately, we will need only the ensemble 
average of the second order velocity, ( 2 Vk). It turns out that in our case of a straight initial field, 
°b a /3 = const, which we consider here, this average is zero, 

( 2 V ka (t)) = 0, (60) 

(Malyshkin 2001). The reason for this simple result is that °b a p is constant in space, and the 
Kolmogorov turbulence, U, is statistically homogeneous. Therefore, the ensemble averages of the 
terms in the brackets [...] in equation (34) are constant in space, and their spatial derivatives are 
zero. As a result, the ensemble averaged velocity ( 2 V^ka) = ( 2 ^ka) is also zero. 



We again assume that u> rot T v ~ 1/3 <C 1. Next we choose such time interval At » tv, that A<f> ~ Lu TO tAt <C 1. 
Then we can simultaneously expand the cos A<j> in eq. (58) and use the random- walk approximation to estimate A(f> 2 . 

16 In particular, note that the driving force (41), which is oc k 2 , becomes small on large scales. 
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To conclude this section, let us integrate equation (54) in time, and obtain the formulas 
rt rt ' T ~- r °° to 2 + q 2 sin 2 (wt/2) du _ J 



[ dt'f H F (t'-t")dt" = ^ f 

J0 J0 K J-i 



UJ 2 + §2 1 + T U} UJ 



'Ok 
2 



t - r(l - e~*/ r ) 



+ 



ft' 



g 2 t - 1 + e" 



-92* 



r 2 t + r 3 (l 



2 g 2 2 



9? 



/•< ft ft ft 

/ dt'/ H F (t'-t")dt" = 2 dt' H F (t' -t")dt" Jofe^t, 



T§*> (61) 



(62) 



which we will use below. Here, the integrals over to can be done by closing the integration contours 
in the complex plane and by evaluating the residues. The final answers in these formulas, written 
after the right arrows "— give the results in the limit t > r, q^ 1 . 



4. Energy Spectrum of Random Magnetic Fields 

In this section we use equations (52), (53) and (60), which give the statistics of the turbulent 
velocities, to derive the evolution of the magnetic energy in the magnetized turbulent dynamo 
theory. 



4.1. The Growth of the Total Magnetic Energy 

The volume averaged and ensemble averaged magnetic energy per unit mass is 

where p is the plasma density, and -B 2 k=o is t ne k = Fourier coefficient of the magnetic field 
strength squared, B 2 . To find £(t), it is convenient to introduce a symmetric tensor B a p = f 
B a B/3 = B 2 b a p. The differential equation for B a @ follows from equation (19), 

dtB a p = B a d t Bp + Bpd t B a = V aa Bp~ f + Vp^B ai - VyB a/3>J . (64) 

Now, we simultaneously solve this equation and equation (20) by making use of the quasi-linear 
expansion procedure, described in Section 3. First, we write 

B 2 (t) = V + X B 2 (t) + 2 B 2 (t), B ap {t) = °B a p + l B a(5 {t). (65) 

Second, we substitute these expansion formulas into equations (20) and (64). We find that the first 
order equations are 

8 t l B 2 = 2 l V a /B a( ]- 1 V ( ?B 2 t(3 , (66) 

d^Bap = l V an °B^ + l Vp tl B ai - ( 1 V 7 °B a/ g) )7 , (67) 
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and the second order equation for 2 B 2 (t) is 



3 t 2 B 2 = 2 l V a /B aP + 2 2 V a /B af3 - fVpWXp - (V/i? 2 ), 



/3 • 



(68) 



Third, we integrate equation (66) in time (with the zero initial conditions) and ensemble average the 
result. Using equation (52), we obviously obtain ( B 2 (t)) = 0. Fourth, we integrate equation (67) 
in time, and then Fourier transform the result in space, r — > k. We have 



(69) 



k' 

k"=k-k' 



Here, we use B a p = B 2 b a p and °b a p = const. (This last formula is our first working hypothesis). 
Fifth, we integrate the second order equation (68) in time. Then we ensemble average the result 
and Fourier transform it in space, setting k to zero. Using equations (53), (60) and (69), we obtain 



'B\ =0 (t) 



2i 



o 



E^O^o^-k^'))^' 



2°B\ =0 I dt' 

Jo jo 



t' 



Ys^^V^it'YV-^t'^dt'' = 2V k=0 ^/i 2 fc 2 



dt' 



t' 



H F (t' - t"; n + ft rd , n rd ) + H F (t' - t"; Q + ft rd , 2U + ft rd ) 



= 2 7 tV 



k=0 



where 



n/n Id = 6k 2 /kl , 2n/n rd = 90 {k 2 /k 2 u ) ^ 2 (i - ^ 2 ) . 



l + (l + 20/O rd )" 2 dfi, 



dt" 
(70) 

(71) 
(72) 



Here, we also use equation (61) in the limit t>r, replace the summation over k by integration, 
making use of d 3 k = 2irk 2 dkd/j,, and use equations (48), (59). 

Now, we derive the differential equation for the averaged magnetic energy £. Following Kulsrud 
and Anderson (1992), we choose t small enough for the quasi- linear expansion to be valid, but large 
enough for the limit t » r to be satisfied. This is very similar to an assumption that the turbulent 
velocities are <5-correlated in time (Kazantsev 1968), the assumption generally used in the dynamo 
theories (see the discussion in the end of this section). As a result, using equations (65) and (70), 
we obtain 

dt (b\ =0 ) = \ (B\ =0 (t) - B\ =0 (0)) = j ('BVo + 2 S\ =0 ) = \ ( 2 S\ =0 ) = 2 7 °5 2 k=0 , (73) 



and using equation (63), we finally obtain 

dS/dt = 2 7 £. 



(74) 
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According to this last equation, the magnetic energy grows exponentially in time, the same 
way as it does in the kinematic dynamo theory (Kulsrud & Anderson 1992). However, the growth 
rate 7, given by equation (71) in the magnetized turbulent dynamo case, is different from the 
growth rate 7 G in the kinematic dynamo case. The latter can easily by obtained by taking the limit 
r2 r d — ► 00 in equation (71) and by changing the viscous cutoff wave number in equation (13) to the 
standard one, k v = (Uo/kov) 3 / 4 ^- In this limit, the back-reaction velocities are zero (because they 
are totally damped), the turbulence is Kolmogorov, and equation (71) reduces to the corresponding 
formula of Kulsrud and Anderson (1992), as one might expect. 

The integrals in equation (71) can be carried out numerically (Malyshkin 2001). The result is 

7 ~ 8.5 {UvL/v) 1 ' 2 (U /L), 7" 1 ~ 10 6 yrs (£/10) 1/2 (L/0.2 Mpc) 3 / 2 , (75) 
7 Collapse ~ l(L/U ) « 10 3 (£/10)- 1/2 (M/10 12 M )" 1/2 , (76) 
7/7o ~ 10. (77) 

Here, i C oiiapse is the protogalaxy collapse time, £ is the ratio of the total mass M to the baryon mass, 
L is the protogalaxy size, and the temperature is assumed to be virial. Equation (77) predicts that 
the magnetic energy growth rate in the magnetized dynamo theory is up to ten time larger than 
that in the kinematic dynamo theory. Two different effects contribute to this difference. First, the 
effective viscosity in the magnetized dynamo case is smaller than the molecular viscosity, v c g = u/5. 
This effect makes the growth rate larger by a factor of y/E (this factor was included in Kulsrud 
et al. 1997). The rest of the contribution comes from the local anisotropy of the turbulent velocities 
in a strongly magnetized plasma. 

Note that according to equations (12) and (75), the magnetic field growth time, 7 _1 = 
0.12(UqL I V) -1 / 2 (L /Uq) , is approximately 40% smaller than the eddy turnover time on the viscous 
scale, t{K) ~ 0.18(U L/v)-V 2 (L/Uo). 17 Thus, the quasi-linear expansion of the MHD equations 
in time may not be fully compatible with the assumption of the <5-time correlation of the turbulent 
velocities. On the other hand, the effects of the finite velocity correlation time should decrease the 
magnetic energy growth rate by a factor of order two (this reduction was found in the kinematic 
turbulent dynamo theory by Schekochihin Sz Kulsrud, 2001). As a result of this reduction, the 
expansion should become better justified. Thus, including the finite time correlation effects into 
our theory would be important but not vital, since our calculation predicts a very large energy 
growth rate anyway [see eqs. (75)-(77)]. 



In the kinematic dynamo theory r(fc„) ~ (1/3)7 1 (Kulsrud & Anderson 1992). 
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4.2. The Mode Coupling Equation for the Magnetic Energy Spectrum 

The ensemble averaged magnetic energy spectrum is 

M(t ' k) = 4^ (^) 3 / k " k)|2) ^ (78) 

where the integration is carried out over all directions of k = k/fc, and Bk is the Fourier coefficient 
of the magnetic field, B. The total magnetic energy, given by equation (63), is obviously 

1 f°° 

£ = - M(t, k) dk. (79) 

2 Jo 

To find the evolution of M(t, k), we use the quasi-linear expansion formula (22) for B. We write 
the ensemble averaged square of the magnetic field Fourier coefficient, up to the second order, 

(|B k (t)| 2 ) = |°B k (t)| 2 + K 1 ^)) ^) + cc] + (rBk(t)| 2 ) + [( 2 B ka (t))°B* ka (t) + c.c] 

= |°B k (t)| 2 + (rBk(t)| 2 ) + [( 2 B ka (t))°BL(t) + cc]. (80) 

Here, c.c. is the complex conjugate, and we use the fact that the ensemble averaged first order 
magnetic field is zero, ( 1 Bj { ) = 0, [this follows from eqs. (27) and (52)]. We calculate the two last 
terms in the right-hand-side of equation (80) separately. 

We start with calculation of the (| 1 Bk(t)| 2 ) term. First, we integrate equation (27) in time, 
and then Fourier transform the result in space, r — > k. We have 



^k x (*) = ikj(6 xa 6 7 s - 5ja8 x s) V] 1 V k " a {t') B k/s dt', 



(81) 



k' 

k"= k-k' 



where we use the divergence free conditions k' a °B k / a = and k'^Vwu = 0. Using this equation 
and its complex conjugate, we obtain 

(| lB k(*)| 2 ) = k-yk T (S xa S jS - S ya S xS )(S x/3 S TV - S T/3 S XV ) 

x f f E E °sw°^» f| < i ^ a (o 1 vtv /J (t ff )>A'A ff 

Jo Jo k , k ,„ 

k"=k-k' k iv =k-k"' 

= k 2 fc 7 fc r (5 aj a(5 7< 5(5 rr; — 5 art 5/3 T 5js — 5 ai 5(i55 Tri + 5 ai 5(3 T 55 ri ) 

x E /7 < Vl Bk'| 2 ( 1 ^k»a(0 1 ^k» /3 (^' , ))^' , • (82) 
,,/ Jo Jo 



k' 

k"=k-k' 



Here, we use ( 1 Uk"a 1 ^k iv /?) <^k",k iv ) see equation (53), and therefore, k lv = k" and k w = k'. 
We also assume that °b a/ 3 = const (our first working hypothesis), and therefore, °B k 's°B k , ri 
k' 



% v \°B kt \ 2 = V|°B-.' 2 
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Now, according to formula (78), we need to integrate equation (82) over all directions of the 
unit vector k. We carry out this integration in Appendix A, the result is 

k 2 {^B^ 2 ) d 2 k = t J dk'IC(k,k') J k' 2 \°B kl \ 2 d 2 k', (83) 



—J J do sin 3 ej ok „j (i + n" /n'l d ) 



" i r,n \ 2 f k ~\- 2k(k - k' cos 9) cos 2 p 



k" 2 



l-(l + 2^/^ d )- 



k" = (k + k - 2kk' cos 



k 2 {k' -k cos ey + (k 2 -k' 2 ) sin 2 9 sin 2 p 2 



k" 2 (k' - k cos 9) 2 + k 2 sin 2 9 sin 2 p> 



,1/2 



ti"/n'; 



VI /rill 

'rd ' 1,2 ' """ ' ""rd 1,2 

ft,,, ft,,, 



cos 2 } , (84) 
(85) 

fc 2 (k' - kcos9) 2 + /c 2 sin 2 sin 2 09 9 „ 9 , . 

90 — ^ sin 2 6» cos 2 <p, (86) 



k" 2 



where function J fc is given by equation (13). 

Next, we calculate the ( 2 B ka (t)}°B ka + c - c - term of equation (80). First, we integrate equa- 
tion (32) in time, and then Fourier transform the result in space, r — ► k. We have 



'B kr) (t) = ik T {5 v p5 TX - 8 Tl3 5 vx ) [ V l Bu x {t') + 2 V k „p(t') °5 k , J dt', 

J o L J 



(87) 



k' 

k"=k-k' 



where we use again the divergence free conditions, k a ®B ka = 0, k a l B ka = 0, fca'v ka = and 
fc a 2 7ka = 0. Second, we ensemble average his equation. The second term in the brackets [...] 
averages out because of equation (60). Third, we multiply the averaged equation by °£? k , add the 
complex conjugate, and use formula (81) for l Bw x . We have 

{ 2 Bk v )°B k + C.C. = ik T (5 v/ 35 TX - Srpdrjx) *(<V^7<5 ~ ^a5 x s) 



£ K f dt ' f dt " fV&ain X V V .p{H)) °Kr?Bu„ 5 + C.C. 

Jo Jo 



k' k'" 

k"=k-k' k iv =k'-k"' 



V (fc 7 - fdt' f dt" ('V-^ait") 'Vvpit')) + C.C. 

Jo Jo 



10 Jo 



= -2fe a %|°Bk| 2 V / dt' I dt"( 1 V^ a (t") 1 V^ f3 (t')) 
k „ Jo Jo 

+ 2k T (6 aT 5 v p - 5 av S TP ) \ r/ \°B k \ 2 V ^ /V /* CV-v^t") 

v , Jo Jo 



r \3 (<oo /-t z-4' 

= - 2fe 2 |°B k | 2 ( —J J d 3 k"j^ dt' J dt" k^V^t") l Vu,p{t'))k p . 
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Here, we use ( ^ k iv a ^ k"/?) oc S k w _ k «, see equation (53), and therefore, k lv = — k" and k"' = k. 
We use the field divergence free condition, k a °Bk a = 0, this is why there are only three terms in 
the third line of equation (88). On the fourth line of equation (88) we use k' = k — k" to change the 
summation over k' to the summation over k". To obtain the fifth and sixth lines of equation (88), 
we use k a °B ka = 0, fc„ 1 V r _k"a = and °B kr? °B k5 = 6 T;<5 |°i? k | 2 (because °b r)5 = const). The 
two terms in the sixth line of equation (88) cancel each other because of the symmetry of tensor 
( 1 V r _ k /' Q , l V\< i iip) with respect to the exchange a <-> j3 [see eq. (53)]. We change the summation over 
k" to the integration over k" in the last line of equation (88) . 

Note that k is perpendicular to °b because the field is divergence free, fc Q °B kQ , = 0. Therefore, 
to carry out the integration over k" in equation (88), we use a spherical system of coordinates, 
<i 3 k" = k" 2 dk" sin 9 d9 dp. Here 6 is the angle between k" and k, so that k" • k = cos 9, and 
ip is the angle between °b and the projection of k" on the plane perpendicular to k, so that 
fi" = k" • °b = sin 6 cos tp. 18 Using equation (53), and equation (61) in the limit t S> r, we obtain 



fdt' f dt" k^V-^ait") l Vu l( ]{t'))kf] 

Jo Jo 



-J ok "t (1 + tt"M' d ) sin 2 e\i- 1 - (1 + 2^7^)' 



cos 2 6 cos 2 ip 



cos z 



' cos 2 <p + sin tp 



Now, we substitute this equation and formula d 3 k." = k" 2 dk" sin 6 d6 dip into equation 
integrate the result over all directions of the unit vector k. We have 



(89) 
and 



k 2 



( 2 B^)°Bi v + c.c. 



d 2 k = -tk 2 



j A; 2 |°B k | 2 d 2 k ( AY fVv, 



V2vr J 



„dk" ( si 



2tt 



x j ckp(i + n"/n" d ) \ i- i - (i + 2n"/n'l d ) 



cos 2 9 cos 2 ip 



Here, Q." ', Q" and Q" d depend on k" and on fx' 



cos 2 9 cos 2 ip + sin 2 ip 
sin 9 cos <p, see equations (48) and (59). 



sin 3 9 d9 

(90) 



Finally, we substitute equations (80), (83) and (90) into equation (78) for the magnetic energy 
spectrum M(t,k). We choose t small enough for the quasi-linear expansion to be valid, so that 
dt(M(t,k)} = [M(t,k) — M(0,k)]/t. As a result, we obtain the mode coupling equation for the 
magnetic energy spectrum, 



dM 



= [ K(k,k')M(t,k')dk' -2^k 2 M(t,k), 
Jo 47r 



(91) 



where the mode coupling kernel K(k, k') is given by equations (84)-(86), and the turbulent diffusion 
constant 

g = \ k"' 2 Jok»i dk'" £ sin 3 9 d9 dip (1 + fV"Kd) 2 



3 Note that these angles 9 and <p have no any relation to the angles in equation (84). 
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«(l-[l-(U20» W 1 2 r 2 ! COS '% }■ 02) 
[ L J cos cos 2 if + sir 99 J 

fi'"/^ = 6^2-, 2fi"7fi^ = 90-^- sin 2 cos 2 99(1 -sin 2 cos 2 <p). (93) 

The function Jq^ is given by equation (13). Here, we use equations (48) and (59). To obtain 
equation (86) we use equations (A8). To obtain equation (93), we use formula p" 2 = sin6>cos</?. 
In the last three equations for the turbulent diffusion constant we replace k" by k'" in order to 
distinguish it from k" in the equations (84) for the coupling kernel. Note, that according to the 
mode coupling equation (91), the k" modes of the turbulence interact with the k' modes of the 
magnetic field to change the energy in the k modes of the magnetic field. 

In the kinematic turbulent dynamo case, we take the limit f2 r d — > 00. In this limit, as one 
might expect, after integrating over <p, equation (84) reduces to the mode coupling equation of 
Kulsrud and Anderson (1992). Equation (92) for the turbulent diffusion constant, after integrating 
over both and (p, reduces to the corresponding equation of Kulsrud and Anderson (1992) [see also 
Kraichnan and Nagarajan (1967), who obtained the same equation in a different form]. 



4.3. The Magnetic Energy Spectrum on Subviscous Scales 

Equation (91), which gives the evolution of the magnetic energy spectrum in the magnetized 
turbulent dynamo theory, is the principal result of this paper. However, this equation is rather 
complicated for easy interpretation. In this section we limit ourselves to the evolution of the mag- 
netic energy spectrum on small subviscous scales, which is of great interest during the magnetized 
dynamo stage in a protogalaxy 19 . In this limit, k 3> k u , and the integro-differential equation (91) 
for the magnetic spectrum simplifies to an ordinary differential equation. 

Let us refer to equation (84) for the mode coupling kernel K(k,k'). Function Jq^" cuts off 
at the viscous wave number k v [see eq. (13)]. Therefore, in the large-fc limit, k 3> k u , we have 
k" ~ \k — k'\ <C k, k', and can expand the kernel K(k, k'). However, the simplest way of calculations 
is to introduce an arbitrary function of k, F(k), which varies slowly in the region k 3> k u , and 
vanishes outside of this region (Kulsrud & Anderson 1992). To derive the mode coupling equation 
on small (subviscous) scales, we calculate the following integral 



F(k) °£dk = ± (±Y r F d -^l *k = J-(L\r F (|B k (*)| 2 )-|B k (0)| 2 ^ 

3 



lirp\2irj J_< 

tp- \ { S^J^ <l lfik WI 2 > d " k + fj (k) [< 2 ^a(i))°£L + c - c -] d3k } 



19 The time evolution of the magnetic energy on large scales is slow because it is set by the turnover time of large 
turbulent eddies. As a result, the large-scale field is not of significant interest until the small-scale field saturates, the 
magnetized dynamo stage ends and the inverse cascade stage begins. 
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F(k) 



,2d 2 M , A 



<9M 



(94) 



Here, to obtain the first line, we use equation (78) and replace the time derivative by the time finite 
difference, assuming that t is small, and our quasi-linear expansion is valid. To obtain the second 
line, we use equation (80). The derivation of the final result [the 3rd line] is given in Appendix B. 
The constants T, Ai and Aq are 



r = 



2\2n J 

r2ir 



k 4 J 0k dk / d9 sin 3 9 cos 2 9 
o Jo 



A, 



o 



-3 + 



dip (i + n/n Id y 
r(^) 



i 



1- (l + 20/O rd ) 



cos 2 9 cos 2 (p 



cos 2 9 cos 2 ip + sin 2 pi J ' 



(95) 



2tt 



fc 4 J k dk / d9 sin 3 9 dp (l + fi/firdj \ 2 

cos 9 cos p + - sin 



i - (i + 2n/n rd )" 



2 + i 



sin 



2 cos 2 9 cos 2 + sin <p 

2tt 



cos 2 cos 2 p 



(96) 



A:' 



/•7T of ] 

J ofc dA; y d9 sin 3 #y dp (l + fi/fird) |2sin 2 0cos 2 <p - - sin 2 9 



+ 



n/n T 



6 



1 - (l + 20/O rd ) 

2r2/i7 r( j = 



cos 2 9 - sin 2 + 



sin 2 9 cos 2 



2 cos 2 cos 2 + sin p> 



cos <£> 



A: 2 



and function J fc is given by equation (13) 



90 -2 sin 2 6> cos 2 p (1 - sin 2 cos 2 99) , 

fC,. 



(97) 
(98) 



Equation (94) is valid for an arbitrary function F. As a result, the integrands in the left- and 
right-hand-side of this equation should be equal, and we finally obtain the mode-coupling equation 
for the magnetic energy spectrum M(t, k) on small (subviscous) scales 



dM 



,d 2 M fk ^BM 

(Ai - l)k— +A M 



Ok 2 v ^ Ok ' " " (99) 

The constant T and the constant dimensionless numbers Ao and Ai can easily be calculated nu 
merically. The result is (Malyshkin 2001) 

r » 11 (U L/u) 1/2 (Uq/L), Ai « 2, A « 5, 



and therefore, 



dM 

~~bT 



,d 2 M 



k 2U1 " k dM 



5M 



(100) 



(101) 



(It is easy to see that Ai = 2 and Ao = 5 correspond to the limit f2 rd <C i^fc 2 .) 

If for the moment we consider the kinematic turbulent dynamo, then we take the limit f2 rd — ► 00 
in equations (95)-(97). After integrating over p and 9, we have F = j , where 7 Q is the Kulsrud- 
Anderson (1992) magnetic energy growth rate, A x = 3 and A = 6. In this case equation (99) 
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coincides with the corresponding equation of Kulsrud and Anderson (1992) and with the corre- 
sponding equation of Kazantsev (1968), as one might expect (note that Kazantsev's equation is 
given in the Fourier space). 

Now, assume that M(t,k Te {) is known as a function of time at some reference wave number 
k = k re {. Then the solution of (99) is 

M(t,k)= I M(t',k iei )G(k/k iel ,t-t')dt', (102) 

J — oo 

where the Green's function G(k, t) is 

G(M ) = ,/jL fcAl/2 ^ e (T/5)(Ao-A;/4)t e -5 In' fc/4Tt = e (4T/6)t g -5 In* fc / 4 rt_ (1Q3) 

V 47r Y L l 2 t 6 l 2 V 47r T y l 2 t 6 l 2 

[This Green function can be obtained by applying the Laplace transformation in time, t — > s, 
to equation (99), see Kulsrud & Anderson 1992; Malyshkin 2001]. To obtain the final result in 
equation (103), we use equations (100) for Ai and Ao- 

According to equations (102) and (103), we see that a "signal" M(t, & re f), at zero time, will 
increase exponentially as e^/^ Tt and will extend down to the scale /c pea k ~ e( 4 / 5 ) r * fe re f , where 
^peak is the peak of function kG(k,t), (of course, the field scale can not become less than the 
resistivity scale). As a result, in the magnetized dynamo theory the magnetic energy tends to 
quickly propagate to very small subviscous scales, the same way as it does in the kinematic dynamo 
theory (Kulsrud & Anderson 1992; Schekochihin, Boldyrev, & Kulsrud 2002). However, the tail 
of the magnetic energy spectrum on fc rc f < k < fc pca k scales increases with the wave number as 
oc fc Al / 2 = k instead of oc k 3 / 2 in the kinematic theory (Kulsrud and Anderson 1992). Note, 
that according to equations (75) and (100), the growth rate of the Green's function, (4/5)T, is 
approximately equal to a half of the growth rate of the total magnetic energy (the latter one is 27) . 



5. Discussion and Conclusions 

In this paper we have developed a theoretical basis for the magnetized turbulent dynamo, which 
operates in protogalaxies, where the plasma is fully ionized, and the viscosity is the Braginskii tensor 
viscosity. The results of the kinematic dynamo theory, already seem to support the primordial 
(protogalactic) dynamo origin of cosmic magnetic fields (Kulsrud et al. 1997). The results that we 
have obtained for the magnetized dynamo, further support this primordial origin theory. This is 
because the number of the magnetic energy e-foldings during the magnetized turbulent dynamo 
stage in a protogalaxy, given by equation (76), is up to ten time larger than that in the kinematic 
dynamo theory. This number of e-foldings is more than large enough for the magnetic field in 
the protogalaxy to grow from its seed value, provided by the Biermann battery, up to the field- 
turbulence energy equipartition value. The number of e-foldings of the magnetic energy on the 
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viscous scale, which is determined by the growth rate (4/5)r of the Green's function (103), is less 
by one half, but it is still sufficiently large 20 . 

Another our prediction is that the tail of the magnetic energy spectrum on small subviscous 
scales increases with the wave number as oc k [see the Green's function (103)], instead of oc A; 3 / 2 
in the kinematic theory 21 . Therefore, in the magnetized dynamo the magnetic energy is slightly 
shifted to larger scales as compared to the kinematic dynamo case. 

The Green's function solution (103) indicates that in the magnetized dynamo theory the mag- 
netic energy tends to quickly propagate to very small subviscous scales, similar to the kinematic 
dynamo case. On the other hand, the observed cosmic fields have rather large correlation lengths. 
Therefore, the magnetic field lines must be unwrapped on small scales by the Lorentz tension forces, 
while the field energy is transferred and amplified on larger scales during the inverse cascade stage. 
This most important and most interesting stage happens when the magnetic field energy is compa- 
rable to the turbulent kinetic energy. In the final part of this paper let us discuss the significance 
of the Braginskii viscosity for the inverse cascade in more details. 

First, note that the theory of the inverse cascade in a plasma with the regular isotropic (non- 
Braginskii) viscosity has a difficulty of unwrapping of the small-scale magnetic field lines. This 
difficulty can be understood as follows 22 . The equation for the turbulent velocities V in the 
plasma with the isotropic viscosity, including the Lorentz forces, is (Landau & Lifshitz 1984) 

d t V a = -P, a + f a + uAV a + (l/4Trp)(B-V)B a -(V-V)V a (104) 

[compare with eq. (4)]. We can estimate the unwrapping velocity, Vunwrap; by Fourier transforming 
equation (104) in space, r — > k, and then balancing the viscous and the inertial forces against the 
magnetic tension force. The isotropic viscosity dominates on small scales, and we have 

^V unwrap ~ (1/47^,1 5 2 , y unwrap ~ (k\\/k) (k u /k) {Vl/uK) « Va, (105) 

where Va is the Alfven speed. At the time of the field-turbulence energy equipartition on the viscous 
scale, the Alfven speed is Va ~ vk v , and before this equipartition Va < vk v . Since in the kinematic 
dynamo theory the field lines have a folding pattern, k S> fen ~ k u (see Fig. 2), the unwrapping 
velocity (105) is small compared to the Alfven speed, even at the equipartition time. In other 
words, since Va ~ V at the energy equipartition, then k\\VA ~ 7 (here 27 is the magnetic energy 
growth rate), and the unwrapping rate, fey Vu nwrap , is much smaller than 7. This means that the 
field continues to grow on the viscous and subviscous scales even beyond the energy equipartition. 



20 Of course, our results (71)-(77) for the magnetic energy growth rate are sensitive to the value of the physical 
parameter fi rd , which is estimated in equation (59). We also left out the finite time correlation effects. Therefore, our 
result for the number of magnetic energy e-foldings should be viewed as an estimate, valid within a factor of order 
two. However, it is important that the number of e-foldings, which we found, is a very large number. 

21 This our prediction is not sensitive to the value of fi rd , as long as fi rd <C vkl. 

22 Cowley, Kulsrud, & Schekochihin 2001, private communications. 
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However, in the case of the magnetized turbulent dynamo the viscosity term in equation (104) is 
modified, and the anti-unwrapping argument does not apply. Indeed, the field unwrapping velocity 
is parallel and varies in the direction perpendicular to the magnetic field lines. The large velocity 
gradient perpendicular to the field lines, which leads to a large perpendicular stress in the isotropic 
viscosity case, is suppressed in the case of the Braginskii viscous forces (because the transfer of the 
ion momentum in the perpendicular direction is inhibited). Therefore, in the magnetized dynamo 
theory fcyVunwrap ~ k\\V\ ~ 7 at the equipartition, and the magnetic field strength saturates on 
the viscous and subviscous scales. As a result, the Braginskii viscosity makes the inverse cascade 
of the magnetic energy more likely, because the larger turbulent eddies do not need to deliver their 
energy to the field on the viscous and subviscous scales (Kulsrud 2000). 
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Goodman, Robert Rosner, Alexander Schekochihin, David Spergel and Samuel Vainshtein for very 
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the Department of Astrophysical Sciences at Princeton University for financial support. This work 
was partially supported by the DOE under the ASCI program at the University of Chicago and 
under DOE Contract No. DE-AC 02-76-CHO-3073. 



A. Details of the derivation of the mode coupling equation 

Equation (82) has four terms in the right-hand-side. Therefore, we have 

^(^Bk^l^^k = T — T' — T" + T'", (Al) 



where there are also four terms 

, 3 coo r r rt rt 



T= ^ 4 (t|) J™ dtf J &\ Bv\ 2 #\{ J ftfkj J C^'ait'^K'andt'dt", (A2) 

T' = k 4 (J^j ' j™ dk 'j k ' 2 l°^k'| W J ^d 2 k^ <o 6 Q ( 1 U k » a (t0 1 Uk» /3 (^ ,/ ))%d^ / d^ /, , (A3) 
T"=T'*=T', (A4) 
T"> = k* (J-J J^dk'J k' 2 |°£ k ,| W| d 2 kjj'* M^k-atO^k-MO)^ dt'dt", (A5) 

Here, k" = k — k', /j, = (k ■ °b) [see eq. (44)], and we replace the summation over k' by the double 
integration over k! and k'. Next, refer to Figure 3A. Vector k' is perpendicular to °b because the 
field is divergence free, &4°^k'a = 0. The following useful equations are valid (Malyshkin 2001), 

k-k' = cos#, k ■ k" = (A; — k' cos 8) / k", k" = (k 2 + k' 2 - 2kk' cos 6) 1 / 2 , (A6) 
fi = k ■ °b = sin 6 cos 92 , <i 2 k = sin 9 d6 dip , (A7) 

« n,- en k „ „ 9 (k' - k cos 9) 2 + k 2 sin 2 8 sin 2 p . 
ti" = °b • k" = — sinflcos^ , 1 - fi" 2 = ± ^ -■ (A8) 
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Using these formulas, equation (53) and equation (62) in the limit t S> r, it is straightforward to 
calculate the double time integral terms in equations (A2)-(A5), 



t r t 



o Jo 



fvv'atfyvwaif)) dt'dt" = j ok »t (i + n"M'd) i + (i + 2n'7nJd) 



-2 



1 + n'7^ d \ 2 k'(k' -k cos 0) sin cos p 



JO JO l fe 



fc //2 

A;' 2 sin 2 i 



(A9) 
, (A10) 



k cos 9) 2 + A; 2 sin #sin 



sin 2 + (1 + 2n'7firi) 



" lr\H \ ~ 2 ~~ ^ COS COs2 V 



fe" 2 



(All) 



Here, 0", fi" and f2" d depend on k" and on fi" 2 [see eqs. (48), (59)]. In turn, k" and //" are 
functions of fe, A;', and [see eqs. (A6), (A8). Now, we substitute formulas (A7) and (A9)-(All) 
into equations (A2)-(A5). The factors, which we obtain in these equations after integration over 
d 2 k = sin# d6 dip, depend only on k and k', so they can be exchanged with the integrations over 
d 2 ~k'. Combining the results together in equation (Al), we finally obtain equations (83)-(86). 
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Fig. 3. — Relative position of vectors °b, k', k" and k = k' + k" in space. Note, that k' is 
perpendicular to °b because the magnetic field is divergence free. The plot on the left (A) refers 
to Appendix A, the plot on the right (B) refers to Appendix B. 
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B. The derivation of equation (94) 



First, we use formula d 3 k = dk d 2 k, equations (A1)-(A5) and equation (88) to obtain the term 
in the brackets {...} in equation (94), 



roc r roo r r 

{...} = J F(k)dk J k 2 (\ 1 B k (t)\ 2 )d 2 k + F{k)dk J k 2 [( 2 B ka (t))°B ka + c.c. 



d 2 k 



7rrl rrll I frill I rj-lV 
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li 2 k 2 F{k)\ Q B k ,\ 2 d 3 k'd 3 k I l { 1 V k » a {t') 1 V k „ a {t"))dt'dt" 

x ../ — x JO JO 

oo roo rt rt 



/oo roo rt rt 

\°B k ,\ 2 d 3 k' ^ 2 k" 2 F(k)d 3 k" Cv^ait'^V^t'^dt'dt", (B2) 
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b a { 1 Vu, a {t') 1 vl,, p {t"))kpdt'dt h 



T" = T' 



T' 



oo J — oo 
oo 



/oo rt rt 

lx"k"F{k)d 3 k" / / °b a ( 1 V k ,i a (t') 1 vt t i (t''))k'pdt' 'dt", (B3) 
-oo Jo JO 

(B4) 

F(k) \ Bv\ 2 d 3 k'd 3 k [ [ k a { l Vu, a {t') 1 Vu, l} {t"))k l 3dt'dt" 
Jo Jo 

/CO roo rt rt 

\°B k ,\ 2 d 3 k' F(k)d 3 k" / k'^V^^t'fV^^k'pdt'df, (B5) 
-oo J — oo JO Jo 

% v = - 2 f°°F(k) \°B k \ 2 d 3 k r d 3 k" fdt' f dt" k a CV_ k//a (t") 'Vv.pWkp 

J — oo J — oo J Jo 

/OO POO f*t ft 

F(k)\°B k \ 2 d 3 k d 3 k'' k a { l V^ k ,, a {t' l ) 1 V kll p{t'))kpdt'dt''. (B6) 

-oo J-oo Jo Jo 

Here, to obtain the final results in equations (B2)-(B5), we change the integration over k in these 
equations to integration over k", using k = k' + k". We use k'^Vwai an d therefore, k a l V k n a = 
k'JVwa- We also use fik = °b • k = °b • k" = k"(°b ■ k") = ^"k" = k" sin 9 cos (p (see Fig. 3B). To 
obtain the final result in equation (B6), we use equations (53) and (62). Next, we calculate T , T£ 
and T"' separately, up to the second order in k" <C k, k' . Refer to Figure 3B. 

First, following Kulsrud and Anderson (1992), in equations (B2)-(B5) we expand F(k) in 
k" <C k at point k' up to the second order. We have, see Figure 3B, 



k = k' + k" cos 9 + (k" 2 /2k') sin 2 9 , 



(B7) 



F(k) = nk') + ^cos9 + ±^k" 2 sin 2 9 + l^k'' 2 cos 2 9, (B8) 

(Kulsrud & Anderson 1992; Malyshkin 2001). 

Second, we calculate %, given by equation (B2). Because fi" 2 k" 2 is of the second order in k", 



-32- 



we need to keep only the zero order term in expansion (B8) for F(k). Thus, we have 

/oo fOO ft ft 

F(k')\°B^\ 2 d 3 k' y" 2 k" 2 d 3 k" / { l Vu, a {t>fK,, a {t"))dt>dt" 
-OC J OO JO JO 

pOO p pOO pit 

= t F(k')dk' k' 2 \°B k ,\ 2 d 2 k' k" A J ok "dk" d9 sin 3 9 
Jo J Jo Jo 

x dip cos 2 <p(i + n"/n'l d ) 2 1 + (1 + 2n"/n't d )~ 2 , (B9) 

Jo L J 

where we use d 3 k" = k" 2 dk" sin 9 d6 dip, n" = sin 9 cos tp and equation (A9). Here and below, 
functions Q", Q" and tt" d depend on k" and y!' 2 = sin 2 9 cos 2 tp, see equations (48) and (59). 

Third, we calculate TJ, given by equation (B3). Because fi"k" is of the first order in k", we 
need to keep only the zero and the first order terms in expansion (B8) for F(k). We have 

/oo fee ft ft 

F(k')\°B k ,\ 2 d 3 k' y"k"d 3 k" / % a { 1 V^ a {t') 1 K„p{t"))k' p dt'dt" 
-oo J -oo Jo Jo 

+ r ^|0^ k ,|2 d 3 k / r^ k ,,2 cos e d 3 k « t f% a ^v^ a {t') i vi, l p{f))k' p dt>dt" 

J -cc.dk J_ 00 Jo Jo 

fOO J Tp r /"OO rn 

= - t k'^j dk' / k' 2 \°B V | 2 d 2 k' / k" 4 J ok , dk" / d6 sin 3 9 cos 2 9 

/•27T 

x / ^cos 2 ^(l + ^7^ / d ) 2 (l + 2^7^ / d )~ 2 . (BIO) 

Here, we use equations (53), (62), k' a °b a = 0, ci 3 k" = k" 2 dk" sin 9 d9 dtp, n" = sin 9 cos 99. The first 
term in the first line of equation (BIO) vanishes after the integration over 9 because the integrand 
is an odd function of cos#. 

Fourth, we calculate T"' , given by equation (B5), keeping all terms in expansion (B8) for F(k), 

coo foe ft ft 
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Here, the first term in the first line is equal to minus T™, given by equation (B6). To calculate the 
other three terms (on the 2nd, 3rd and 4th lines), we again use equations (53), (62), k' a °b a = 0, 
<i 3 k" = k" 2 dk" sin 9 d9 dip, /j," = sin 9 cos (p. The term in the second line of equation (Bll) vanishes 
after the integration over 9 because the integrand is an odd function of cos6>. 

Next, we substitute equations (B4), (B9), (BIO) and (Bll) into equation (Bl). The — term 
of equation (Bll) cancels the T" term in equation (Bl). Then, we substitute the result, obtained 
in equation (Bl), into the second line of equation (94) and, using equation (78), we find 
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(B14) 
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Now, we integrate the right-hand-side of equation (B12) by parts over some extent in k! and choose 
F(k'), so that it and its derivative dF/dk' vanish at the end points. We finally obtain the third 
line of equation (94), by dropping the double primes, ", in equations (B13)-(B15), and introducing 
the new constants T, Ai and Aq, given by equations (95)-(97). 
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